Comment on "Equation of state of dense and magnetized fermion system" 
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Contrary to what is claimed by Ferrer et al. [Phys. Rev. C 82, 065802 (2010)], the magnetic field of a neutron 
star cannot exceed 10^'' G and the thermodynamic pressure of dense magnetized fermion gas is isotropic. 
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The authors of recent paper |1] construct thermodynamics 
of charged fermions in strong magnetic field B where Lan- 
dau quantization of orbital motion is important and thermody- 
namic quantities depend on B. The subject attracts consider- 
able attention, with the most important applications to neutron 
stars possessing strong magnetic fields. The authors conclude 
that (i) the magnetic field in a neutron star can exceed 10^^ G 
and (ii) the gas of particles in a quantizing magnetic field has 
anisotropic pressure. We point out that both statements are 
inaccurate. 



I. MAXIMUM FIELD STRENGTH 

The well known estimate based on the virial theorem |2] 
gives the maximum magnetic field in a neutron star Bmax ^ 
10^® G |i3i] . The authors of Ref. {V\ claim that this estimate can 
be relaxed. As an alternative, they propose arbitrary simplis- 
tic parametrizations of mass density p and field strength B as 
functions of the radial coordinate r within the star, treating the 
parameters of these functions as "totally arbitrary." For certain 
values of these parameters they obtain Smax > 10^^ G. 

However, the density and field distributions are not arbi- 
trary, but must satisfy stability equations for a magnetized star 
with a realistic equation of state. Detailed self-consistent nu- 
merical simulations (for example, |4, 5]) show that, depending 
on the adopted equation of state in the stellar core, Bmax takes 
values (0.3 - 3.0) x 10^^ G, in disagreement with Ref. [IJ but 
in good agreement with Ref. \%. A large variety of equa- 
tions of state were explored in numerical simulations j^]. The 
obtained p and B distributions are different from artificial dis- 
tributions of Ref. Ol, leading to different values of Bn-^a^. 



II. ISOTROPY OF PRESSURE 

The consideration of the pressure in Ref. Ol] is based on 
the articles by Canuto and Chiu who showed that kinetic 
pressures p|™ and p*^'" of charged particles along and across 
B, calculated as ensemble averages of respective currents of 
kinetic momenta, are different. The authors of Ref. lUl] repeat 
the consideration (§] using a more general formalism and ar- 
rive at the same conclusions. According to Refs. 011^, the 
total anisotropic pressure is the sum of the magnetic pres- 
sure related to the Maxwell stress tensor, and the kinetic pres- 



sure. The longitudinal and transverse kinetic pressures are 
p'^'" = —fl ~ MB and py'" — -il, where fl is the grand 
canonical potential per unit volume and M is the magnetiza- 
tion (directed along B in the quasistationary approximation 
adopted in these studies). 

However, the deficiency of the approach of Ref. Oal has 
been pointed out long ago by Blandford and Hernquist ^ . It 
is well known that the total microscopic electric current den- 
sity j is composed of the free (or conduction) current term 
jf and bound current term jb due to magnetization (dynami- 
cal polarization contribution to jf, in the quasi-stationary ap- 
proximation is negligible). The magnetization cuiTent den- 
sity equals (in Gaussian units) jf, = cV x M; in case of 
boundaries, this volume current should be supplemented by 
the surface current cM x B/B (see, e.g., Ref. |8]). The 
total thermodynamic pressure P in a magnetized plasma is 
the sum of the kinetic pressure and an additional contribution 
due to the Lorentz force density related to the magnetization 
cuiTents. If we compress a plasma across B, then the mag- 
netization cuiTent density induces an additional contribution 
MB to the force density. As a result, the transverse com- 
ponent of the total (thermodynamic) plasma pressure equals 
p'^'" + MB = Py'", so that the total plasma pressure P = — fl 
is isotropic. 

In spite of simplicity of the above arguments, they are 
sometimes ignored in the literature, like in Ref. 1 1]. Therefore, 
in order to make them still more transparent, let us illustrate 
the pressure isotropy with two graphic examples. 

As the simplest example, consider a plasma contained in 
a finite cylinder in vacuum with a uniform external B-field 
along the cylinder axis. At equilibrium in the absence of ex- 
ternal forces, the sum of the force densities exerted on the 
side wall of the cylinder by transfer of kinetic momenta of 
plasma particles and by the surface magnetization current 
equals p*^'" + MB = —fl. It is the same as the force den- 
sity p|™ = — exerted on the head wall. Hence the plasma 
pressure, which can be determined in this experiment by mea- 
suring forces on the cylinder walls, is isotropic. 

As another example, more relevant to astrophysics, con- 
sider a volume element in a magnetized star. Let the ele- 
ment be sufficiently small and distributions of B, tempera- 
ture T, and gravitational acceleration g be sufficiently smooth, 
so that we can assume constant J3, T, and g within this 
volume. Let the z axis be directed along g. Then p and 
il{p,B, T) depend on z, resulting in z-dependent magnetiza- 
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tion TVf = —dVl{p,B^T)/dB. Hydrostatic balance implies 
the density of gravitational force, pg, be balanced by the den- 
sity of forces created by plasma particles (gradient of kinetic 
pressure and Lorentz force due to plasma magnetization). 

Now let us compare two limiting cases. If B is parallel 
to g, the z-component of Lorentz force is absent, and we 
get the standard equation of hydrostatic equilibrium pg = 
= dP/dz = -dfl/dz. 

If B is perpendicular to g, then the kinetic pressure gra- 
dient dp^ /dz acts in parallel with the Lorentz force density 
BdM/dz. Note that in our case dM/dz ^ 0, simply because 
dp/dz ^ (p depends on z) in the gravity field. Since B and 
T are constant. 



dM _ dM{p,T,B) dp 
dz dp dz 



d^n{p,T,B) dp 



dz 



dp dz dpdB 
Then the equilibrium condition takes the same standard form 

dnif'" dM d , ^ , dM dfl 

pg=JL^+ B— MB) + B— = - — . 

dz dz dz dz dz 

Thus, the gradient dp/dz = —{dil/dp)^^pg does not de- 
pend on B-field direction, which means that the hydrostatic 
equilibrium is determined by the isotropic thermodynamic 
pressure P, in accordance with the results of Ref. |01. 

Since the forces created by bound currents are small in the 
majority of appUcations, the equations of magnetohydrody- 



namics (MHD) are commonly derived neglecting the magne- 
tization. However, the magnetization term is easily recovered 
by substituting the general expression j ~ jf + jb into the 
microscopic Lorentz force density j x B/c that is included 
in the derivation of MHD equations from the first principles 
(e.g., 10], Chap. VIII). Moreover, thermodynamics of magne- 
tized media is well studied in the theory of magnetics (e.g., 
|9], Chap. IV). Of course, everyone is free to use anisotropic 
kinetic pressure in MHD equations and add the magnetization 
force density explicitly. However, it seems more natural to 
follow the traditional approach and use the isotropic thermo- 
dynamic pressure that automatically includes the contribution 
of the magnetization. 
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